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1. Introduction 

In this paper we rigorously analyze the validity of a numerical technique for studying 
resonances in quantum mechanics. The technique is called, "A spherical box approach 
to resonances," by its inventors, Maier et al ||. We prove that the technique yields 
correct energies and lifetimes for shape resonances in the semiclassical limit. 

The technique is an "L 2 method," in contrast to time-independent scattering theory 
methods, such as the calculation of phase shifts near energies where a resonance is 
expected. These L 2 methods are surveyed, e.g. , in ||. 

The basic physical idea underlying all L 2 methods is that a resonance wavefunction 
is a state that is concentrated mainly in the interaction region. In contrast, states 
associated with the rest of the continuous spectrum are not concentrated in any bounded 
interval. As a consequence, when the system is confined to a box that is large compared 
to the interaction region and the size of the box is varied, the resonance wavefunction 
is much less influenced than the states from the rest of the continuous spectrum. This 
should be visible in the spectrum, and is the basis of the technique we study. 

To make this precise, we consider the Schrodinger operator 

H:=D 2 + V, D:=-^- (1) 

i ax 

with a resonance producing potential V that is defined on all of M. We restrict the 
system to the interval (—£, £) with Dirichlet boundary conditions at x = ±£, and plot 
the eigenvalues of the resulting operator H (£) as a function of I. 

Figure [I] presents the results obtained by doing this for the potential V that is 
depicted in Figure 0. 



E(lj 




I 



Figure 1. An example of the dependence of the eigenvalues on the box size £. 
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In this figure, one can clearly distinguish between eigenvalues that depend strongly 
on £ and others that seem to be almost independent of £. Furthermore, there are avoided 
eigenvalue crossings when a strongly dependent eigenvalue is close to an eigenvalue that 
is almost independent of £. Note that in our example, eigenvalues are not expected to 
ITJ , since the potential has no apparent symmetry properties. 



cross 



In addition to relating the almost constant eigenvalues to resonance energies, 
Maier et al f| also relate the sizes of the gaps in the avoided crossings to the imaginary 
part (or width, or inverse lifetime) of the resonance. In |§, spherically symmetric 
potentials are treated. After the reduction to an angular momentum subspace, the 
particle can escape to infinity in only one way, by increasing the radial coordinate r. 
In the model we consider, the particle can escape toward either plus or minus infinity. 
Since the probabilities for going in the two directions can be different, we observe two 
different size gaps for each given resonance. This is obviously the case in Figure |l|. For 
our model, the resonance width is related to the larger of the two gaps. 

In this paper, we provide rigorous justification of these results in the semiclassical 
limit. As a first step, we adopt a standard definition of a resonance that is presented 
Kj . This definition identifies a resonance with a complex eigenvalue of a suitably 



m 



constructed analytic family of operators obtained from the original Hamiltonian (p]). 

In many instances, as in the case of shape resonances, such a complex eigenvalue 
can be viewed as arising from the perturbation of an eigenvalue embedded in the 
continuous spectrum. We take this viewpoint and employ the framework of, "The 
Shape Resonance," M by Combes et al. We temporarily impose supplementary Dirichlet 
boundary conditions at points u± to decouple the interaction region from the rest of 
M. This yields an unperturbed operator on all of M that has embedded eigenvalues 
whose eigenfunctions are supported in the interaction region. Removal of these Dirichlet 
conditions perturbs the embedded eigenvalues to produce the resonances (that are 
realized as complex eigenvalues of certain non-self-adjoint operators). The perturbation 
calculations are facilitated by the use of Krein's formula [||. 

To relate the resonances of H defined on L 2 (M) to the almost ^-independent 
eigenvalues of H(£), we show that the techniques of j|] can also be applied in a box to 
study H(£). We then employ the following strategy: For small values of h, resonances of 
H are very close to embedded eigenvalues of H with supplementary Dirichlet conditions 
at u±. For £ > max{|u; + |, |^-|}, these embedded eigenvalues are also eigenvalues of 
H (£) with supplementary Dirichlet conditions at uj±. For large £ and small h, removal of 
these supplementary Dirichlet conditions perturbs these eigenvalues only slightly. Thus, 
the resonances of H are near eigenvalues of H(£). These results are made precise in 
Theorem |^. 

This approach also allows us to prove rigorously that the gap in the avoided crossing 
is on the order of the square root of the resonance width, in accordance with ||. We 
accomplish this by relating both the gap and the resonance width to the thickness of 
the potential barrier as measured by the Agmon distance |T[]. The relationship between 
resonance widths and Agmon distances is already established in [|], so we need only 
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examine the relationship between the Agmon distances and the gaps in the avoided 
crossings. This is done in Theorem |3|. 

2. Hypothesis and Results 

For simplicity, we assume the potential V to be bounded. We wish to study resonances 
that are produced by a single well and to avoid asymptotically degenerate eigenvalues 
with an exponentially small separation in h. Furthermore, we want the bottom vq of 
the well to be above the scattering threshold. We force this situation by imposing a 
hypothesis that can be expressed nicely with the help of the notion of the classical 
forbidden region at energy E. This is defined as 

J(E) :={xeR: V(x) > E} . 

Our precise hypothesis is the following: 

(HI) V G C 1 (-K) is bounded and has a local minimum vq at Xo, such that J(vq) is 
connected, and limsupui^ V(x) < v . 

By translating the origin if necessary, we choose an interior region 

Qi := (u;_, u)+), with c<j_ < and uj + > 0, such that fli \ {x } C J(v ). 

We define the exterior region to be Q e := M \ Q i} and let Vt~ = (— oo, ujJ) and 
= (u) + , oo). We define the decoupled comparison operator H d as having the same 
symbol as H, but with supplementary Dirichlet conditions at and u> + . This operator 
decomposes into 

H d = H l ® H e with V(H a ) = HlnH 2 (n a ), where a G {i, e} . 

Since we want to focus on shape resonances, we impose a hypothesis that prevents 
resonances from being produced in the exterior region for energies near vo- We phrase 
this hypothesis in terms of a non-trapping condition || : We say the potential V is non- 
trapping in Q e at energy E (abbreviated E is NT), if the following condition is satisfied 
for a G { — , +}: 

3S > 0, Vx G O" \ J(E), ^—^(2(V(x)-E)+xV'(x))<-S. (2) 

We assume: 
(H2) v is NT. 

Note that formula (H) is equivalent to the more standard virial condition 



3S > 0, VxeO e \ J(E), 2(V(x) - E) + xV'(x) < -S . 
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Figure 2. The potential associated to Figure 1 and relevant parameters 

Furthermore (|2|) implies the "exterior" virial condition 

3S > 0, Vi6^ \ J(E), (2(V(x) - E) + (x - u a )V'(x)) < -S, a G {-, +} . 
Our third hypothesis concerns analyticity under exterior dilation. For 9 G M, we define 

U e :L 2 (R)-^L 2 (R) by 

{uo_ + e e (x — lu__), x < lo^ 
x, x G (u;_, cu + ). 

cj + + e e (a; — cj + ), x > uj + 

We then assume: 

(H3) Vq := UqVUq 1 defined initially for 9 E Ft, has an analytic continuation as a bounded 
operator to the strip {9 G(F : |Im#| < /3}, for some (3 G (0, 7r/4). 

For 9 G -K we also define the operators #0 := UeHU e l and ifg := UgH d U e 1 . It is a 
straightforward calculation to obtain the associated symbol 

Ug(D 2 + V)^ 1 =r'- 2 D 2 + Vor , 

—h 2 u"(x), x G (id-, uo + ) 



where [r' e 2 D 2 u](x 



-fi 2 e~ 2d u"(x), x £ w + ] . 



Since t/g is a unitary operator on L 2 (M) for 9 <E M, we easily compute the domains for 
the operators Hq and Hq, for 9 G iR: 



Z>(#e ) = D(#*) © £>(iT 



V{H g ) = {m ®u e e H 2 {tti) © H 2 {Q e ) : u e (w±) =e5Ui(w±), u' e (w±) =e^u' i {u ± )}. (3) 
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We define the restrictions of these operators to the box (—£, i) to be 

H${£) ■= r'~ 2 D 2 + Vor e on V{Hf) n ((-■£, *)), and 

:=r'i- 2 D 2 + Vor e on V(H e ) n J)) . (4) 

For = 0, Hq=o(£) is simply the Schrodinger operator if (£) described in the introduction 
that is used to produce plots, such as Figure 1. 

The following lemma describes the analytic continuations of these families of 
operators to complex values of 6: 

Lemma 1 Hypotheses (HI) (H3) imply the following two conclusions: 

(i) {Hg, \lm8\ < (3} and {Hg(£), \lm8\ < (3} are self-adjoint analytic families of Type 
(A) of m- sectorial operators. 

(ii) {Hg, \lm9\ < /3} and {Hg(t), \lm9\ < (3} are self-adjoint analytic families of 
operators. 

Proof: These conclusions for the families Hg and Hg are proved in Q. The same 
proofs apply for the families Hg(£) and Hg(£) since the proofs in j| make no use of the 
(un)boundedness of fl e . a 

We next recall the Agmon distance Q, that we denote by the symbol ds- It is the 
distance associated to the pseudo-metric ds 2 := max{0, V(x) — E}dx 2 . We introduce 
the abbreviations 

d v '■= d v (xo, oc£), a G {-, +}, and d* := min{^ , cf+} . 

The following theorem gives precise information about the resonance on the line 
and the "resonance in the box." Its first conclusion follows from 0]. 

Theorem 2 Assume (HI) (H3) and that E d is the n th eigenvalue of H % . 

(i) For any $ G (0, 1) and sufficiently small h~, there exists (3q G (0,(3), such that H^ 
has a ( complex) eigenvalue E close to E d that satisfies: 

E = E d + y^— ^, with t = o(e- 2M * /h ), where a n = o(l), Vn > 1. 
' n! 

n>l 

(ii) The same is true for the operator i7 i/ 3 (£). Furthermore its eigenvalue is stable in 
the sense of Kato /Q, Sec. VIII. 1. 4], as the box size I tends to infinity. As i tends 
to infinity, this eigenvalue converges to the corresponding eigenvalue of H^. 

(Hi) For sufficiently small h and those values of I, for which there exist positive constants 
c and N, such that dist(E d , a(H d (£)) \ {E d }) > ch N , there exists a real eigenvalue 
of H(£) close to E d that satisfies the same type of expansion as above. 
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Remark: (a) Note that in M, the theorem is stated with d* replaced by d Vo (xo, dQi). 
Due to the possible choices of lu±, the difference between the two quantities can be made 
arbitrarily small and can be absorbed into But then, how small % must be chosen, 
depends on 

(b) In the third conclusion of this theorem, one cannot expect uniform results in £ 
and h. The eigenvalues of the exterior operator H e (£) have different dependence on £ 
and h than the eigenvalues of the interior operator H l . The condition dist(.E d , o~(H d (£))\ 
{E d }) > ch N is technical; we cannot handle exponentially closely spaced eigenvalues. It 
is well known that under our hypotheses, the eigenvalues of H l near the bottom of the 
well (close to vq) cannot be spaced more closely than 0(/i 7 ). Here, the constant 7 is 
strictly smaller than 2. Its value depends on how flat the bottom of the well is. In order 
to prove that eigenvalues from H e (£) do not cause dist(E d ,a(H d (£))\{E d }) > ch N to be 
violated for all £, we would need an additional assumption on the decay of the potential. 
For example, together with dilation analyticity, it would be enough to assume that V 
tends to a limit at infinity like \x\~ e for any e > 0. 

We now turn our attention to the gaps in the avoided crossings that occur in graphs of 
the eigenvalues of H(£). For this part of our analysis, we replace hypotheses (H2) and 
(H3) by the following: 

(H4) V G C 3 (M). For i£Sl e \ J(v ), the potential obeys V(x) < v and there exist two 
constants v± < Vq, such that V — v± = 0(\x\~ e ) as x tends to ±00. Furthermore 
for n = 1, 2, we have = (9(|x| _<E_1 ) as x tends to ±00. 

This hypothesis allows us to use WKB estimates to analyze the behavior of eigenvalues 
of H e (£). We note that H e (£) decomposes into the direct sum of H e _{Jl) and H^(£), 
where H e _(£) acts on L 2 ((-£, and H e + {£) acts on L 2 ((u + , £)). 
We have the following result on the gaps: 

Theorem 3 Assume (HI) and (H4)- Suppose E d is an eigenvalue of H l and of 
H^(£ ), but not of Ht a (£ ). Assume it satisfies dist(E d ,a(Ht a (£o))) > ch N , for some 
positive constants c and N and a G { — , +}. Then we have the following: For fixed 
values of h that are sufficiently small, there exists a neighborhood V(£o) of £ , such that 
for all £ in V(£o), H(£) has two eigenvalues E + and E_ that are exponentially close to 
E d . These two eigenvalues are separated by a gap that satisfies 



mm 

tev(^o) 



{\E + -E_\}= 



(h + t 2 ) n a, 



n 



, where a n = o(l), Vn > 1. 



In this estimate, t\ and t 2 satisfy the following for any $ G (0, 1): 



t 1 =t 2 = o(exp(-&d-/h)) if E d ea(H l )na(H e _(£ )), 
tl =t 2 = o(exp(-Mt /h)) if E d ea(H i )na(H e + (£ )). 
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Remark: (a) Here, H does not depend on 

(b) The width of the resonance is given by the tunneling parameter t according to 
Theorem 0. We do not know whether the resonance is going to tunnel to the left or right, 
so we only obtain the estimate t = o(exp(—i)2d*/h)). In Theorem ^ we know to which 
side the resonance escapes, and the result is more precise. We also obtain estimates 
for both of the gap sizes that can occur in the avoided crossings for a given resonance. 
Since d* = min{<i~ , d+ Q }, Theorem || shows that the larger gap is of the same order as 
the square root of the resonance width. We again note that in ||, a radial symmetric 
situation is studied, so that there is only one way for the resonance to escape, and hence 
only one gap size. 

(c) The eigenvalues of H l are obviously independent of £, but not of u±. Thus, it 
might seem that the condition of having a double eigenvalue is crucially dependent on 
the choice of u±. This is not case: From Theorem ^ (iii) we see that the eigenvalues 
of H l vary only by an exponentially small quantity in h when the u± are varied. For 
the eigenvalues of H±(£), we show in |Appendix B| , that (H4) implies that eigenvalues 



E e a(H^(£)) that belong to an interval (t> , v + 6) are related to h, £, and a quantum 
number m by the asymptotic formula 

E = v a + f(m + h^>j 2 (1 + 0(h) + 0{£^)) , a e {-,+}. 

We thus have the following consequence: Suppose, for example, that the n-th eigenvalue 
E d of H l coincides with an eigenvalue of H+(£q) for some choice of u±, and that E d is 
at least a distance of 0(h N ) from the spectrum of Ht(£o)- Then for any other choice of 
u±, there exists an £ in a neighborhood of £q, such that E d is an eigenvalue of H^,(£), 
and the distance from E d to the spectrum of Ht(£) is still at least 0(h N ). 

3. The Proofs 

Inspection of the proofs of for Theorem |2| (i) shows that they are valid whether or 
not Q e is bounded. Furthermore, these proofs can be separated into two parts: The 
first is a formal algebraic part that shows the stability of the eigenvalue of H l for the 
whole operator and constructs the asymptotic expansion of the perturbed eigenvalue in 
powers of the tunneling parameter t. It is quite simple and short. The second part is the 
justification of these algebraic formulas with the corresponding estimates. This part is 
more complicated and involves estimation of the operators involved in Krein's formula. 
We present the formal algebraic part, which is needed in all of the situations treated 



in Theorems Q (ii), (iii), and Theorem |3|. We do this in Section [TT] in the context 
of Theorem |2| (iii). In Section we treat the stability of the resonance eigenvalue of 
Hipo W as the box size £ tends to infinity. Finally, in Section |3.3| we prove Theorem [| 



In the Appendix, we recall Krein's formula and present the more technical estimates, 
including the WKB estimates. 

We omit the estimates required to prove the existence and the series expansion of 
the eigenvalue of Hi/3 (£) because they are identical to those in [§]. 
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3.1. Stability and Tunneling Expansion for The Box 

We view H(£) as a perturbation of H d (£). This perturbation involves two Dirichlet 
conditions. It is most easily approached by way of Krein's formula, that exhibits the 
difference of the resolvents of H(£) and H d (£) as a rank two operator. 

The norm of this rank two operator is not small. However, because the Dirichlet 
conditions are imposed inside the classically forbidden region, its norm does not explode 
in proportion to the inverse of the distance from the spectrum to the spectral parameter 
in the resolvents. This allows us to choose the parameters in such a way that the 
resolvent of the resolvent of H d (£) is small in norm, and we can still use perturbation 
theory. 

The tunneling expansion is based upon a Feshbach type reduction of the eigenvalue 
equation with respect to the unperturbed eigenprojection. This leads to an implicit 
equation that we solve by using the Lagrange inversion formula. 

3.1.1. Stability To simplify the notation, we suppress the £ dependence in many 
of the formulas. We define 

R d (z) := (H d (£) - z)- 1 and R{z) := (H(£) - z)' 1 . 

We choose a contour T that lies in the resolvent set of H d (£) and encloses only E d in 
o-{H d {£)). We then choose a point zq in the intersection of the resolvent sets of H(£) 
and H d (£), but outside of T. By using the identity 

R d (z ) - - J— ) 1 = -(z- z ) -(z- z ) 2 R d (z), (5) 
z- Z J 

we obtain the following expression for the eigenprojection P d = P d {£) associated to E d : 

^ = -ds/ r * w * = -^i(* ( *>- i )~ 1 < s '. 



where {T := z e(T : z = z G T }. By defining 

n(z ) := (H{£) - z^- 1 - {H d {£) - z )-\ 

we can formally write the eigenprojection P = P(£) associated to the perturbed 
eigenvalue E as 

p = I (^o) - ^r 1 (i+^o) {R^-zyy 1 dz. 



If we can choose Y and zq, such that 



< 1, then the inverse term 



n(z )(R d (z )-z) 1 

in the integral for P can be computed by geometric series. Then the eigenprojection is 
well defined, and by standard arguments, we can deduce the stability of the eigenvalue 
for H(£). 
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To see that we can do this, fix any n e IV. Let 

A := dist(£ d , a(H d (£)) \ {E d }), and fix r E [mm{h n , ^A}, -A] . (6) 

Note that by hypothesis, A > ch , for some N e IV, and that we can choose r to be 
as small as any power of H. 

We define r := {2 6(T : \z — E d \ = r} and z = E d + 2ir. Then, formula (|5]) implies 
(i? d (zo) ~z) 1 = 0(r). Thus, the stability follows from the following proposition that 
we prove in Appendix A.l| : 



Proposition 4 n(z ) = 0(1). 

3.1.2. Tunneling expansion Since we have proven the stability of the eigenvalue 
and constructed P(£), we can write the eigenvalue equation as 

R(z )P(£) = -j^-P(i) ■ 
E - z 

We perform a Feshbach type reduction to this equation, with respect to the projections 
P d and Q d = 1 — P d . We define the "reduced" resolvent 

R(z; zo) ■= Q d (Q d (R(z ) - z)Q d ) 1 Q d . 

It satisfies the following estimate: 

Proposition 5 For any z in the disc delimited by T, one has R(j^; Zo) = 0(r). 

Proof : If we replace the R(zq) by R d (zo) in the definition of R(z; Zq), we obtain a trivial 
result. The conclusion to the proposition is obtained by applying perturbation theory 
to this trivial result. a 

For (E — zq)^ 1 the reduction yields the implicit equation 

- Eh~J pdp = pd - «*> *<^v ^ »<*>) pdp - 

Using properties of the trace and the factorization n(z) = hA*(z)B(z), c.f. [Appendix 



A.l| , we obtain 

1 1 .... / ,-. ( . ^ I 



E-z E d -z 
or equivalently 

1 1 



E-z E d -z E-z 
where (suppressing z in A and B) 



frTr (B(z )P d A*(zo)(l - hB(z )R(^- Q , *oM*(*b))) 



ts(- ), (7) 



t := h \Tr(BP d A i ')\ and s(z) := ~ Tr (hBP d A*{l - hBR(z; z )A*)^j . (8) 
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For any z in the disc delimited by T and z = -f— , we have the following estimate on 
s(~z): 

\s(z)\ < \\l-hBR(z;z )A*\\ = 1 + 0(r). 



This follows from Proposition |5] and the bound on tt, c.f. |Appendix A.l| . If we can 



establish the estimate t = o(e 2l9d * / h ) of Theorem then equation (0) can be solved 
with Lagrange's inversion formula || p. 250] 



1 ^ t r ' 



E-z E d -z 

n>l 



d 



71-1 



1 , 1 ^t T - 



E d ~z ' E d -z 4^ n! 



n>l 



Multiplying by {E—Zq){E—Zq) and rearranging, we obtain 

# = ^ - (^-.o)(^-,o) E ^ = ^ d - E(*°-^) fe+1 (E^ 

n>l ' fe>l ^n>l ' ' 

We estimate the coefficients a n by using the Cauchy formula 



We define cr n := Zo)(E d — z )a n and easily obtain the estimate cr n = o(l) of Theorem Q. 

3.1.3. The tunneling parameter The above calculation relies on the estimate 
t = o{e' 2 ^/ h ). To prove this, we note that if <pd denotes the eigenfunction associated 
to E d , then using the definitions and estimations of [Appendix A.l| , 



t < h\\B<f> d \\ \\A<j> d \\ < h 2 \\TRT*\\ \\B<j> d \\ 2 

= " r ^" 2 ^ ^d^-)| 2 + I^K)| 2 )- 

For each part of Theorem |2|, we can estimate the expression \<p' d {u)J)\ 2 + |<^(u; + )| 2 by 
the well known decay estimates of Agmon PJ. This implies the results of Theorem |2|. 

3.2. Stability as The Box Size Tends to Infinity 
We consider the operator 

H?(i):=H e (£)eH?(i), (9) 

where H$(£) is the operator defined in (f|), and 

E e e e [l) := e~ 2e D 2 + Vor on Hi fl H 2 (M \ [-£, £}) . 

It is easy to see that Hq e {£) is an analytic family of Type (A) in 9, and that we have 
the following resolvent estimate: 
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Proposition 6 Assume (HI )-(H3) and let S denote the constant in the non-trapping 
condition (H2). Let v = {z E<£ \ \Rez — u | < S/4, Im z > —S/4}. Then 



Vzev, \\R%(z)\\ < ^(1 + 0(0)). 



Proof : HfS(£) decomposes into a direct sum of operators that act on L 2 ((— oo, —£)) 
and L 2 ((£, oo)). We consider only the term associated to the interval (£, oo); analogous 
formulas hold for the other term. We mimic arguments of ||. For u G 1~L\ D 7i 2 ((£, oo)) 
and any v G L 2 ((£, oo)), we have 

\\v\\\\(H%(£)-z)u\\ > Re((H%(£)-z)u,v). 

For (3 > we use this with v = —ie~ l2/3 u to obtain 

Re {{H^{£) - z)u, v) = - Im [e 2l(: \e- 2l(l 'D 2 + Vor lf3 - z)u, u) 

= - Im (e 2if3 (V o r i/3 - z)u, u) 

= - (((3 (2(V - Rez) + (x-u+)V -Imz) + 0((3 2 ))u, u) 

> {(3 (S - 2(v -Re z) + Imz) + 0{f3 2 )) \\u\\ 2 

S , /n /n2\ \ II 112 



> [^— + 0(P 2 )J \\u 
For negative f3 we repeat this calculation with v = ie~ l2f3 u. This proves the proposition. 



We now fix 9 = i/3 as in [Qj. With the definitions of z and T as in Section |3TX| , we 
define 

PvkW = ~i~ I {(PiM-^r 1 ® {(h^-zo)- 1 -t}- 1 dz. 

Here, Pip {£) projects onto the eigenspace for the eigenvalue E G o-(H i/3o (£)), but does 
so in the space L 2 (M). To prove stability of the eigenvalue in the generalized sense 
(c.f. Kato, ]7], Sec. VIII. 1.4]), we must show that Pip (£) — —> Pip as I tends to oo, 
where 



p i(] = -7T~- {( H ifo-Zo) 1 ~z) 1 dz. 
im Jf 



It is shown in [|J that for sufficiently small h, {{H^ q — Zq) 1 — z) 1 = 0(r), uniformly on 



T. The estimates of [§| are also valid for ((Hip Q (£) — Zq) 1 — z) \ So, from Proposition || 
and identity we see that, 

uniformly on T. Thus, we need only show that for any u G L 2 (M), 
lim || ((H iPo - zo)- 1 - (H° o (£) - zo)- 1 ) u\\ = , 



uniformly in h. This is shown in |Appendix A. 2 
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3. 3. Proof of Theorem 3 

In the degenerate case, we must solve for two eigenvalues. So, we cannot a priori use 
the Lagrange inversion formula to solve equation (0) in the disc delimited by V . 

However, we could use the formula if one of the solutions were known to be — • 

' E d —zo 

This would happen if tt were a rank one operator. In that case, the spectra of H d and H 
would intertwine, and as a consequence, at the crossing of two eigenvalues of H d there 
would have to be an eigenvalue of H. 

In our situation such a scenario can be realized by lifting the two Dirichlet conditions 
one after the other. 

It suffices to consider the case where E d G o~(H % ) fl a(H^_(£ )). In the first step, we 
consider the operators 

H d _(£) := H e _{£) © H* and H-{1) := D 2 + V on L 2 {{-£,u + )) . 

By hypothesis, % is small and fixed, and H l has the eigenvalue E d , which for £ = £q 
is a distance of 0{Ti N ) from the rest of the spectrum of H d (£o), i.e. E d is a simple, 
conveniently isolated eigenvalue of H d (£ ). Thus, the analog Theorem |2| (iii) is valid: 

Lemma 7 Assume the hypotheses of Theorem @ with E d G a(H { ) n a(H%{t Q )). Then 
there exists a neighborhood of V(£o), of size ch N , such that for each £ G V{£q), the 
operator H-(i) has an eigenvalue £L close to E d that satisfies the following for any 
d e (0,1) 

E_ = E d + \" with t = ( e -2*M-->-o)A) and Gn = n\ y n > i 

^— ' n! 

n>l 

Proof : We first note that as we vary £, with the restriction that \£ — £ \ < ch N , E d 
remains isolated from the rest of the spectrum by a distance of size ch . Thus, we can 
prove the lemma by mimicking the proof of Theorem ^| (iii). a 

For the second step, we note that due to the behavior of E+(£) there exists an £% G V(£o), 
such that 

We now use the intertwining of the spectra of H^{£\) © H+(£i) and H{£\). We obtain 
the following lemma by using the techniques we used for Lemma [7] and noting that the 
eigenfunction (f) d associated to E_ has the form 4> d = 0_ ©0+, where if_(£i)0_ = E^(p_ 
and ^(4)0+ = E-<f>+: 

Lemma 8 Assume H(l) and (H4) and that E_ is a double eigenvalue of H_(£i) © 
H+(£i) as constructed above. Then the operator H{£\) has two eigenvalues EL and E + 
that satisfy 

E + = E_ + y ^ tl+t2 ) n(Jn with a n = o(l), Vn > 1 , 
/ — ' n! 

n>l 
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where, for any i9 G (0, 1), 

t x = ( e -2*^(xo^+)/R) anrf £2 = ( e -2*M«+,<i)/R) . 

The last step in the proof of Theorem |3] is to note that the first two steps can 
be done for any admissible u+. The n-th eigenvalue E d of H l changes by only an 
exponentially small amount in h when uj + is varied, so it remains properly isolated from 
a(Ht(£o))- Furthermore, by the behavior of the exterior eigenvalues, there exists an £ 2 
in a neighborhood of Eo, such that the new EL is also an eigenvalue of Hl^)- The 
optimal estimate is obtained when t\ = t 2 , in which case we have t\ = t 2 — o{e~' ddv o/ h ). 

Appendix A. Krein's Formula 

Since we need Krein's formula for one and two supplementary Dirichlet boundary 
conditions, taken at different points depending on the situation, we wish to present 
the formula in a general setting. On the other hand, for simplicity, we leave out the 
exterior dilation. We deal with this only when necessary. 

Suppose n > 2, and — oo < xo < xi < . . . < x n < oo are specified. Let Q := (xq, x n ) 
and Q k := (x k -i,x k ) for k = 1, . . . ,n. Let H := D 2 + V be a Schrodinger operator on 
Q, with self-adjoint boundary conditions at xq and x n , and let H d be the corresponding 
decoupled operator with supplementary Dirichlet conditions at X\, x 2 , ■ ■ ■ , x n -i- Denote 
their resolvents by R and R d , respectively. 

Let z G p{H) n p(H d ) and u,v G L 2 (M). Define u := R d (z)u and v := R{z)*v. 
Clearly, u G V(H d ), and thus, u = Q) k=1 Uk with u(x k ) — 0, k — 1, . . . , n — 1. We have 

((R{z) -R d (z))u,v) = {u,v)-(u,v) 

n 

= ^2(D 2 u k ,v) nk - (u k ,D 2 v) nk 

k=l 

n 

= —h 2 > u k v 

k=l 
k=l 

We use standard Sobolev space notation and define functionals by the following 
relations, where / G ® k=1 'H 1 (Vlk): 

11 rWiilO ■ Tij:= lim f(y), for j = k, k + 1, fc = l,...,n-l. 

If / = T^ k +1 f for all /, we simply write T Xfc . It is well known that is compact, and 
consequently, (TJ )* : (F 2 — > W _1 (fij) is continuous. Furthermore, Lemma 4 of Section 
III of shows that whenever \ G Cq°(M) satisfies x( x fc) = 1 f° r = 1 3 • • • , n — 1, 

||T^m|| 2 < 2h~ l \\xuj\\ \\DxUj\\ < c/T 1 \\xuj\\ n i , for j = k, k + 1. (A.l) 
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( T l \ 






T~ : = 


fey 


, T + : = 




, T d := -T- ©T+ and T := 





With these definitions, we have the following formula, 

((R(z) - R d {z))u, v) = h (R{z)T*iT d DR d {z)u, v) , 
where all the multiplications are understood to be matrix multiplications. 

Appendix A.l. Applying Krein's Formula for Theorem [| (Hi) 

In the proof of Theorem |2| (iii), we have xq = —£, X\ = — u)_, x 2 = £3 = £, 

R d {z) = {H d {£) - z)-\ and R(z) = (H(£) - z)' 1 . 
Following 0|, we define 

B(z) := iT d DR d (z) and A(z) := TR(z) . 
Since H is self-adjoint, we can write 

tt(z) = R(z) - R d (z) = %A*{z)B(z). 
Furthermore, since TR d (z) = 0, we have 

Tn(z) = TR(z) = A(z) = hTR(z)T*B(z) . 
We combine the two formulas to obtain 

tt(z) = Ti 2 B*(z) TR(z)T* B(z) . 
Proposition f| now follows from 

Proposition 9 Let zq = E d + 2ir. Fix any N e IV. Then for sufficiently small h and 
any r £ [min{h N , |A}, |A], 

B(z ) = 0(h 



-l/2\ 



and TR(z )T* = O^ 1 ) 



Proof: The assertion on TR(zo)T* is proved in step 5 of the proof of Theorem III. 3 of 
0. As for B(zq) we have 

ll^o)|| 2 = \\TLDR e (z )\\ 2 + \\T*_DR\z )\\ 2 + \\T? DR(z )\\ 2 + \\T*_DR%z Q )\\\ 



where R^zq) := (W - ^q)" 1 and R e (z ) := (H e (£) - Zo)' 1 . Let X be a C °° function 
supported around u± such that X {oj±) — 1- Using the estimate ( |A.1|) , it suffices to find 
a uniform bound on the expressions 

X DR\z ), D X DR\z ), and X DR e (z ), D X DR e (z ) . 

We choose X , such that V(x) — v > e > for x G supp X . Then steps 1 and 2 of the 
proof of Theorem III. 3 of H show that all theses terms are uniformly bounded. _ 
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Appendix A. 2. Applying Krein's Formula for Theorem |] (ii) 

Here we consider the difference of the resolvents of the operators Hip defined by 
formula (f|) of the introduction and H^ o (£) defined by formula @. In this case, n = 3, 
Xq = — oo, x\ = —£, x 2 = £, and x 3 = oo. The difference of the resolvents is 

R iPo (z ) - Rf^zo) = hR lPo {z»)T*ie- 2 ^T d DR% o {z Q ) . 

Let x be a Cg° function supported around ±£, with x(=t^) = 1- To show that 
T d DR®p Q (z ) and TR_ il3 (z ) are uniformly bounded operators we use the estimate 
( |A.1| ). Thus, it suffices to show that xDR^ (z ) and xR-ip(zo) are uniformly bounded 
operators from L 2 to H, 1 . If that is true, then for u € L 2 (M), we have 

\\(Rip (z ) - Rf^zo)) u\\ < \\TR^(z )\\ (\u(-£)\ + \u(£)\) 

<c(\u(-£)\ + \u(£)\)^0, 

uniformly in h, since u = xDRf^ o (z Q )u G 7i 1 (iR). 

We now address the required uniform bounds. Commuting x an d D, we need only 
consider xDR^ o (z ) and xD 2 R^ (z ). The expressions for TR_ if 3(z ) are analogous and 
can be treated the same way. The formula 

llx^g (^o)|| 2 = ReR° o (z y X 2 {D 2 - 2h 2 { X 2 )") i&(«>) 

shows that it is sufficient to bound xD 2 R®p (z Q ) and xRip ( z o)- We have 

X D 2 R° (z ) = e 2 ^x (H° -z -(Vo r iPo - z )) R° Q (z ) 
= e 2 ^ {1 - (V o n , - z ) X RF, (zo)) 

For xRip ( z o), we set rc if3Q (z ) := (H^^-z)- 1 - {Hfp Q {£)-z)~ l , and then write 
X i2gb(*b) = X (rf(*b) © - ^o)- 1 + © J2&(zb) 

The right hand side is uniformly bounded in h and £ by Propostion|6| and Lemma II. 3 and 
Theorem III. 3 of ||, which are also valid for (Hfo {£) — z ) _1 and iTip (zo), respectively. 



Appendix B. WKB Estimates 

For these estimates, we follow Olver |9|, Ch. 11]. The goal is to find approximate 
solutions to the differential equation 

-h V + (V - E)w = (B.l) 

in Q e with vq + 5 > E > vq for some positive 5. Due to either the non-trapping 
condition or the explicit assumption (H4), there is at most one turning point in each 
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of the intervals (c<j + , oo) and (—00, ujJ). There is exactly one, if 8 is sufficiently small. 
We assume 5 has been chosen so that this is the case. 

It suffices to consider the interval (u + , 00), and we denote the turning point by x t . 
We define a new independent variable £ := s(x) by 

s(x) s'(xf = E - V{x) , s{x t ) = 0, s'{x t )>0. 

By integration, we obtain 

(g \ 2/3 /•max{x,xt} 

-S(x)\ where S(x) := / vT^XO — E\dt . 

" J Jmin{x,xt} 

Note that sgn (V(x) — E) = sgn (x t — x). It is easy to check that under our conditions, 
Theorem 3.1 of ||, Ch. 11] shows that equation ( [B.l|) has two C 2 solutions w\ and w 2 
in (u + , 00), such that 



(B.2) 

Higher order approximations are also known, c.f. |9|, Sec. 11.7]. 

The Dirichlet boundary conditions imply the quantization condition 

Wi(u + ; h)w2{£] H) — w 2 (uj + ; h)wi(£; fi) = . 

Factoring the error in ( B.2|) in the classically forbidden region, using the asymptotic 
expansions of the Airy functions || p. 392/3], and substituting all this into the 
quantization condition yields 



e v " \ cos( 

n 



\) + <>(*)) + i e ^)/» (sin(^ -1) + = . 



If this equation is satisfied, then necessarily, cos(^^ ~ f ) = 0(h). This implies 
^ - f = ^vr + 0(H), or equivalents 

/ ^W^V{fydt= (n + -)Kh + o(h 2 ) . 

J x t 4 

Now using (H4), we have 

/•£ /•£ 

/ ^£ - y(t) = / ^ - v + rft + / (v 7 ^- v(t) - ^E- v+ )dt 

J x t J x t J x t 

V +l * ; Jxt ^E - V(t) + 

= e^/E^(i + o(r e )) 

From this, it follows that 

E=v++ (^±m^ (1+0(h)+0(r , )y 
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